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Abstract 

The ab oscillations in absorbing matter are considered. The standard model based on 
optical potential does not describe the total ab transition probability as well as the chan- 
nel corresponding to absorption of the 6-particle. We calculate directly the off-diagonal 
matrix element in the framework of field-theoretical approach. Contrary to one-particle 
model, the final state absorption does not tend to suppress the channels mentioned above 
or, similarly, calculation with hermitian Hamiltonian leads to increase the corresponding 
values. The model reproduces all the results on the particle oscillations, however it is ori- 
ented to the description of the above-mentioned channels. Also we touch on the problem 
of infrared singularities. The approach under study is infrared-free. 
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1 Introduction 



The theory of ab osciUations is based on one-particle model [1-4]. The interaction of particles a 
and b with the matter is described by a potentials Ua,b- Imf/b is responsible for loss of 6-particle 
intensity. The wave functions '^a,b are given by equations of motion. The index of refraction, 
the forward scattering amplitude /(O) and potential are related to each other, so later on the 
standard approach is referred to as potential model. 

Description of absorption by means of Irallh is at least imperfect and should be partially 
revised. As an illustration, let us consider the case of strong 6-particle absorption. Instead of 
periodic process, we get two-step system decay 



Here {b — medium) — > / represents the 6-particle absorption. The potential model does not 
describe this process as well as the total ab transition probability (see [5,6] and next section). 
It describes the probability of finding a 6-particle only. 

By way of specific example we consider the nn transitions in medium followed by annihilation 



where M are the annihilation mesons. The qualitative process picture is as follows. The free- 
space nn transition comes from the exchange of Higgs bosons with the mass mn > 10^ GeV 
[8,9]. Prom the dynamical point of view this is a momentary process: Tc ~ i/mn < 10~^^ s. 
The antineutron annihilates in a time ~ l/F, where F is the annihilation width of n in the 
medium. We deal with two-step process with the characteristic time T2 ~ Tq. 

The potential model describes the probability of finding an antineutron only, whereas a 
main contribution gives the process (2) because n annihilates in a time r^. In the following we 
consider the nn transitions since the final state absorption in this case is extremely strong. 

So the model should describe the two-step processes like (1) and (2). On the other hand, 
as absorption Hamiltonian tends to zero, the well-known results on particle oscillations should 
be reproduced. This program is realized below. 

Also we present here an elaborated derivation of the lower limit on the nn oscillation time 
and discuss in some detail the uncertainties connected with medium corrections. 

The paper is organized as follows. In the next section, we recall the main results and point 
out the chief drawback of the potential model. In this model the nn transition probability 
depends critically on the antineutron self-energy E. In the field-theoretical approach the similar 
picture takes place (Sects. 3 and 4). Because of this, the particular attention is given to the 
suppression mechanism and origin of S. In Sect. 4 we arrive at the conclusion that S = 0. An 



(a — medium) — > (6 — medium) — >• /. 



(1) 



[7-13] 



{n — medium) — > (n — medium) — > M, 



(2) 
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important effect of tlie competition between scattering and anniliilation of n in tlie intermediate 
state is studied as well. The main calculations are performed in Sects. 5 and 6. If E 0, the 
^'-matrix amplitudes contain the infrared singularity. For solving the problem the approach 
with finite time interval (FTA) [14] is used. It is infrared free. First of all we verify the FTA 
by the example of exactly solvable potential model. The FTA reproduces all the results on 
the particle oscillations {ua'^b: etc). In Sect. 6 the process (2) and process shown in Fig. 
8b are calculated. The linkage between the ^'-matrix theory and FTA is studied as well. In 
Sect. 7 we complete the calculation of process (2). Also we arrive at the conclusion that for 
the processes with zero momentum transfer the problem should be formulated on the finite 
time interval. The results are summarized and discussed in Sect. 8. The limiting cases and 
effects of absorption and coherent forward scattering are discussed as well. Section 9 contains 
the conclusion. 

2 Potential model 

In this section we touch briefly on the main results and the range of applicability of potential 
model in the case of nn transitions. The chief drawback of this model is given as well. 

Let Un =const. and Un =const. be the neutron potential and optical potential of n, 
respectively. The background fleld Un is included in the unperturbed Hamiltonian Hq. The 
interaction Hamiltonian has the form 

Ti-I = 'Hnn + 'H, 
nnn = e^n'^n + H.C, 

n ^ {ReV + iluiV)^n'^n, 
V^Un-Un^ ReUn - Un - iV/2. (3) 

Here Tinn and Ti are the Hamiltonians of nn transition [11] and n-medium interaction, respec- 
tively; e a small parameter e = l/r^jj, where Tnn is the free-space nn oscillation time. 

The model can be realized by means of diagram technique [5,12] or equations of motion 
[10,11,13]: 

{idt — Ho)n{x) = en{x), 
{idt — Hq — V)n{x) — en{x), 

Ho^-Vy2m + Un, (4) 

n(0, x) = 0. For V = const, in the lowest order in e the probability of finding an n in a time t 
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is found to be 



W^{t) = -^[1 - 2cos(Rel^i)e-r*/2 + 6"^*]. (5) 

V r 



In the following we focus on the most important case Ft » 1. Then Wn{t) ~ e^/ | V 1. 
The total nn transition probability Wl°^ (more precisely, the probability of finding an n or 
annihilation products) is given by 

The index "pot" signifies that the non-hermitian Hamiltonian (3) is used. 

The free-space nn transition probability Wf is Wf — eH"^. Comparing with Eq. (6), one 
obtains the suppression factor it!pot: 

wr^ r 1 

The energy gap SE = V leads to very strong suppression of nn transition in the medium 
and changes the functional structure of the result: Wf ~ — Wf"^ ~ t. Because of this, the 
particular attention is given to the suppression mechanism and origin of 6E. 

The energy gap is the antineutron self-energy: 6E = V = S. Indeed, the result (6) can be 
obtained by means of diagram technique: 

Wr' = 2lmTiit, (8) 

= -eG-P^V, (9) 
1 1 



^pot 



en-PrJ2m-Un V 
Here G^°^ is the antineutron propagator, p = (e„, p„) is the neutron 4-momentum; e„ = 
p2/2m + Un. 

Consider now the range of applicability of model (4). The total nn transition probability 
Wtis 

Wt^Wn + Wa, (10) 

where Wa is the probability of finding the annihilation mesons {i.e. the process (2) probability). 
The potential model describes correctly the Wn. However, for the calculation of Wt and Wa 
it is inapplicable. In the one-particle model the total process probability W^"^ is obtained by 
means of Eq. (8) which follows from the unitarity condition. Since the Hamiltonian (3) is 
non-hermitian (in the first approximation one can put ReV — 0), the (S'-matrix is non- unitary 
and Eq. (8) is invahd [5]. The condition of probabihty conservation 

m 
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can be used only if the ^-matrix is unitary or unitarized. 

It must be emphasized that in the problem under study the unitarity of ^'-matrix is of 
particular importance because ImV^ enters the leading diagram and plays the crucial rule (see 
(9)). 

As a result, in the potential model the effect of final state absorption acts in the opposite 
(wrong) direction [6], which is not surprising, since the unitarity condition and non- unitary 
^'-matrix are mutually incompatible. The condition SS'^ = 1 is applied to the essentially non- 
unitary ^'-matrix. (The potential model does not describe the Wt at all. The non-unitarity 
is only formal manifestation of this fact.) As shown below, the potential model also does not 
describe the competition between scattering and annihilation of n in the intermediate state and 
time-dependence of process (2). The greater the |Imy|, the greater an error in the Wf°^ and 
Wa- So Eq. (6) is incorrect. The direct calculation of the antineutron absorption (process (2)) 
is called for. 

3 Free-space process 

First of all we consider the free-space nN annihilation (see Fig. la) and free-space process 

n + N ^n + N ^ M (11) 

shown in Fig. lb. 




Figure 1: (a) Free-space nN annihilation, (b) Free-space process n + N^n + N^M. 
The matrix element of ^-matrix and amplitude corresponding to Fig. la are defined 

as 

iTa =<M| rexp(-i j dxHuNix)) -l\nN>^ Ar„(27r) V(p/ - p^^a- (12) 
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Here < M \ represents the annihilation mesons, 'HnN is the Hamiltonian of nN interaction, A^^ 
inchidcs the normahzation factors of the wave functions. 

Ta and Ma involve all the nN interactions followed by annihilation including the nN rescat- 
tering in the initial state. The same is true for the subprocess of nN annihilation involved in 
Fig. lb: the block should contain all the nN interactions followed by annihilation. 

We write the formulas corresponding to Fig. lb. The interaction Hamiltonian is 

'Hi = Tinn + 'HfiN- (13) 

Formally, in the lowest order in Hnn the amplitude of process (11) is given by 

Here Go is the antineutron propagator. Since Ma contains all the nN interactions followed by 
annihilation. Go is bare. We emphasize this fact as it gives an insight into origin of E. 

Due to the zero momentum transfer in the nn-transition vertex the 4-momenta of n and n 
are equal. The both pre- and post-nn conversion spatial wave functions of the system coincide: 
I -^^p>sp=| >sp- Actually this is true for any neutron state (for any nuclear model). 

For the time being we do not go into singularity Go ~ 1/0. It results from the zero 
momentum transfer in the vertex corresponding to Ti-nn- The value of E is disconnected with 
Tinn and we want to separate these problems. The general consideration is given in Sect. 6. 



4 nn transitions in the medium 

In this section the origin of E (energy gap) is studied in the framework of microscopic theory. 
It is shown that the value of E is uniquely determined by the definition of the annihilation 
amplitude of n in the medium. It turns out that for a realistic definition E = 0. Also we 
consider the competition between scattering and annihilation of n in the intermediate state. 



4.1 S'-matrix approach 

Let us consider the process (2). (The nn transitions with n in the final state arc considered in 
the next section.) We use the scheme identical to that for process (11) with the substitution 
nN — > (n — medium). The background field f/„ is included in the neutron wave function 
(Hamiltonian Hq); the quadratic terms Tinn are included in the Tij: 

Til — Tinn + Ti^ 

Hit) = [ d^xHix). (15) 
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Here H is the hermitian Hamiltonian of n-medium interaction. The sole physical distinction 
with the model (4) is in the Hamiltonian 7i. Recall that the potential model does not describe 
the processes (2) and (11) [5,6]. 

The process amplitude M2 is uniquely determined by the Hamiltonian (15) 

^ 6„ - pl/2m - ^^^^ 

(see Fig. 2a). The matrix element of ^'-matrix Tj?j and amplitude of antineutron annihilation 
in the medium are 

iTf^ =<f\T^ |0n,>= N{2n)'6\pf-p,)M:^, 

/oo 
dtH{t)) - 1 (17) 
-00 

(compare with Eq. (12)). Here | On^ > is the state of the medium containing the n with the 
4-momentum p, <f\ denotes the annihilation products, A'^ includes the normalization factors. 

/ / 

- - " n V - 

\ ! \ 



Figure 2: (a) nn transition in the medium followed by annihilation. The antineutron annihi- 
lation is shown by a circle, (b) Same as (a) but the antineutron propagator is dressed (see 
text) . 

The definition of the annihilation amplitude through Eqs. (17) is natural. If the number 
of particles of medium is equal to unity, Eq. (17) goes into (12). The annihilation width F is 
expressed through M^: F ~ J \ p. Since H appears only in the M^, the antineutron 
propagator is bare. In the next section we perform the rigorous calculation of M2. 

It is important that M2 ~ M^. The value of F and corrections to (if they are possible) 
have little effect on the results. 



7 



Constnict now the model with the dressed propagator (see Fig. 2b). In the Hamiltonian H. 
we separate out the real potential V — Ret^ — 



(18) 



and include it in the antineutron Green function 



G 



G"^ + G'^VG^ + ... 



1 



1 



(19) 



(l/G^) - V 



V 



Then 



M2 = eG'^M' 



(20) 



(21) 



The propagator G"* is dressed: H — V 0. According to (21), the expressions for the 
propagator and vertex function are uniquely connected (if Hr is fixed). The "amplitude" 
M'^{V,H') should describe the annihilation. However, below is shown M'^ and model (20) are 
unphysical. 

We recall the amplitude Ma involves all the nN interactions followed by annihilation in- 
cluding rescattering in the initial state. Similarly, involves all the n-medium interactions 
followed by annihilation including the antineutron rescattering in the initial state. Compare 
now the left- and right-hand sides of (21). 

Prom the physical point of view model (20) has no justification on the following reasons. 

1) If the number of particles of medium n is equal to unity, the model (20) does not describe 
the free-space process (11) because Eq. (14) contains the bare propagator. 

2) The observable values (F for example) are expressed through M™ and not M^. Compared 
to M™, is truncated because the portion of the Hamiltonian H is included in G"^. has 
not a physical meaning. 

(The formal expression for the dressed propagator should contain the annihilation loops as 
well. In this case the statements given in pp. 1) and 2) are only enhanced.) 

3) Equations (19) and (20) mean that the annihilation is turned on upon forming of the 
self-energy part T, = V (after multiple rescattering of n). This is counter-intuitive since at the 
low energies [15-17] 



and inverse picture takes place: in the first act of n-medium interaction the annihilation occurs. 

The realistic competition between scattering and annihilation should be taken into account. 
Both scattering and annihilation vertices should occur on equal terms in M™ or G™. According 
to pp. 1) and 2) the latest possibility should be excluded. In line with the physical meaning of 
and M2, the amplitude (16) allows for the above-mentioned effect (see Sect. 4.3). 



r = 




> 2.5 



(22) 
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The structure with dressed propagator hke (20) arises naturally if V and Ti.' are the prin- 
cipally different terms and vertex function does not depend on V. In the problem under 
consideration this is not the case. This is evident from the formal expansion of the T-operator 

Texp(-z J dx{V^n'^n + n')). (23) 

It is significant that even the non-realistic model (20) gives reinforcement in comparison 
with the potential model [5]: 

Wa , , fr/2 



>1 (24) 



because for the model (20) the probability of process (2) is 

W^a ~ r (25) 

instead oiWt" ~ l/F. 

To summarize, the introducing of dressed propagator G"^ (energy gap) into process model 
entails an uncertainty of the vertex function M^. The all-important effect of the competition 
is not taken into account. The limiting case n = 1 is not reproduced. is unknown and 
unphysical. 

We do not sec the reasons for existence of field V which should be included into G"" and thus 
the antincutron propagator is bare. For the process sown in Fig. 8b the propagator is bare as 
well. Essentially, this fact is governing. Below we assume the definition (17) and consequently 
the model with the bare propagator (16). 



4.2 Simplest model 

The fact that the antineutron propagator is bare is obvious in the model containing the anni- 
hilation vertex only. We consider Fig. la. Assume that 

= ^M^a^n^AT, (26) 

where $m denotes the fields of mesons. The diagrams of the nN annihilation are shown in Fig. 
3. 

Similarly, for the n-medium annihilation we take 

^ = K = E*M^«^-^^.- (27) 

i 

The corresponding diagrams are shown in Fig. 4. 

Consider now the process (2) using the same Hamiltonian Tia- The diagram is shown in 
Fig. 5; the Hamiltonian is given by Eqs. (15), where 7i = 7ia- The antineutron propagator is 
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Figure 3: nN annihilation. The interaction Hamiltonian is given by Eq. (26). 




Figure 4: Antineutron annihilation in the medium. 

bare. The questions connected with the self-energy part do not arise in principle because Ti^ 
must appear only in the TJ^ (see Fig. 4). The block TJ^ is described by Eqs. (17) and (27). 

In view of Eq. (22), the models like (27) are reasonable and so it seems obvious that the 
antineutron propagator is bare. 

4.3 Scattering and annihilation of n in the intermediate state 

In the low-density limit the relative annihilation probability of the intermediate antineutron ri 
is [15-17] 

ri = — > 0.7, (28) 

= CTa + Cs, where a a and Ug are the cross sections of free-space nN annihilation and nN scat- 
tering, respectively. The ratio (28) or (22) is very important for the correct model construction. 
The model given above reproduces the magnitudes of r and r-i. Indeed, let us consider the 
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N 




Figure 5: nn transition in the medium followed by annihilation, 
free-space process 

n + N ^n + N ^ f, (29) 
where / denotes nN or M. The annihilation and scattering channels are defined by (11) and 

n + N ^ n + N ^ n + N, (30) 

respectively. The corresponding diagrams are shown in Figs, lb and 6a. Using the amplitude 
(14), the cross section of process (11) is found to be 

= N J d<^\ Mib |2= a^N J d^\Ma P= aV„, (31) 

a — eGq. The normalization multiplier N is the same for o"^^ and a a- 




" x" 




" x" 








N 




N 



Figure 6: Free-space processes n + N^n + N^n + N (Fig. a) and n + N^n + N^n + N 
(Fig. b). 

For process (30) the similar calculation gives 



nN 2 



a (To 



(32) 
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and correspondingly 



^ = - = ^- (33) 



The model (13) reproduces the ratio r. 

For nn transition in the medium ri will be calculated by means of optical theorem. To 
check this calculation we obtain ri for free-space process (29) by means of optical theorem as 
well. The on-diagonal matrix element (see Fig. 6b) is 

M(0) = eG'oM,(0)G'oe = a'M,(0), (34) 

where Ms({)) is the zero angle nN scattering amplitude. Let o""^ be the total cross section of 
process (29). Using the optical theorem in the left- and right-hand sides of (34), we get 

a,"^ = aV* (35) 

and 

^ = - = (36) 

For process (29) the relative probability of the annihilation channel is given by (28), as we 
wished prove. 

In the medium instead of (11) and (29) one should consider the process (2) and inclusive 
nn transition 

(n — medium) ^ (n — medium) — /m, (37) 

respectively. Here fm denotes M or n. The result is the same (see Appendix A): for process 
(37) the relative annihilation probability of the intermediate n is given by (28). 

Ratio (28) is exphcitly used only in the classical models hke cascade one [18]. However, the 
Hamiltonian should contain all the needed information, which allows the calculation of r or 
ri. The fact that the model reproduces these ratios is very important; otherwise one can get 
a wrong, additional suppression as in model (20). Since the potential model does not describe 
the processes (11) and (2), it cannot reproduce (33) and (79). 

The principal results of Sect. 4 are as follows, (a) The antineutron propagator is bare 
and singular, (b) In the low-density limit the ratio (28) should be reproduced. This can be 
considered as a necessary condition for the correct model construction. Model (15) satisfies 
this requirement. 



5 Field-theoretical approach with finite time interval 

The model must satisfy the following requirements: a) The S'-matrix should be unitary, b) 
The model should reproduce the free-space process shown in Fig. 1 and competition between 
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scattering and annihilation considered above. These conditions are obvious, however they are 
not fulfilled in the potential model. 

The interaction Hamiltonian is given by Eq. (15). We use the basis (n, n). The results 
do not depend on the basis. A main part of existing calculations have been done in n — n 
representation. The physics of the problem is in the Hamiltonian. The transition to the basis 
of stationary states is a formal step. It is possible only in the case of the potential model 
H — V =const., when the Hamiltonian of n-medium interaction is replaced by the effective 
mass H rUes — ReV — iV /2. Since the calculation of process (2) will be done beyond the 
potential model, the procedure of diagonalization of mass matrix is unrelated to our problem. 

The yS-matrix amphtudes corresponding to Figs, lb and 2a are singular as Go ~ 1/0 and 
~ 1/0. Contrary to quantum electrodynamics, the formal sum of series in e gives the 
meaningless self-energy E ~ e^/0. This is because the Hamiltonian T-tnn corresponds to 2-tail. 
There is no compensation mechanism by radiative corrections. 

For solving the problem the FTA is used [14]. It is infrared free. The calculation is performed 
by means of the evolution operator U {t, 0). The limiting transition i — > oo is not made as it is 
physically incorrect. The FTA can be used for any problem since for the nonsingular diagrams 
it converts to the .S-matrix approach (see Sect. 6.1). 

5.1 nn transitions with n in the final state 

First of all we consider the nn transitions with n in the final states on the finite time interval 
(t,0) (see Fig. 7). 
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Figure 7: nn transition in the medium with n in the final state. 
We introduce the evolution operator U{t,0) — I + iT{t,0). In the lowest order in e the 
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matrix element is given by 

<nO| U{t,0) -I\On>= iT^n(t,0) = -i <npO\ [' dtcHnn{tc) + T^{t,0) dtcHnniQ \ Onp>, 

Jo Jo 

(38) 

ft °° ft ftk-i 

T^'it^Q = rexp(-z / dt,H{t,)) - 1 = ^(-z)' / dt,... / dtkH{ti)...H{tk), (39) 

Jtc k=l 

where | Orip > and | Orip > are the states of the medium containing the neutron and antineutron 
with the 4-momentum p— (e^, Pn), respectively; T" is the T-operator of n-medium interaction 
(compare with (17)). 

We expand the 'I'-opcrators in the cigenf unctions of unperturbed Hamiltonian Hq = — V^/2m+ 
Un- Taking into account that Hnn \ On^ >= e | Ofip >, we change the order of integration [14] 
and obtain 

Tnn {t,0) = -et-e [' dtciT^i {t-Q, (40) 
Jo 

iTt,{T) =<ripO|T^(T) |Ori^>, 

where r — t — T^{t) = T^\t,tc). The n-medium interaction is separated out in the block 
T^{r). This equation is important since the structure of matrix element corresponding to the 
process (2) is similar (see (64)). On the other hand, Eq. (40) can be verified with the use of 
the exactly solvable potential model. 



5.2 Verification of FTA 

To verify the FTA we obtain the results (5) and (6) of the potential model. As in Sect. 2, we 
take H — V =const. The block T^{t) is easily evaluated, resulting in 

iTt,{r) = C/^(t) - 1 = exp(-iyT) - 1. (41) 

The probability of finding an n is 

Wn{t,0)=\T^n{t,0)\'' . (42) 

By means of Eqs. (40) and (41) it is easy to verify that | Tjj„(i, 0) p coincides with Eq. (5). 
The total nn transition probability Wt is given by 

Wt{t, 0) - 1- I Unit, 0) 21mTii{t, 0), (43) 

where Uii{t, 0) =< UpO \ U{t, 0) | Onp >. In the framework of the FTA the on-diagonal matrix 
element Tu has been calculated in Ref. [14]: 

Tii{t, 0) = ieHy2 - f dta dUTJl{ta - Q. (44) 

Jo Jo 
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Using Eqs. (41) and (44), one obtains that 2lmTu = Wf°*. 

Consequently, the FTA reproduces all the results of the potential model. This was to be 
expected since one and the same Hamiltonian was used. The same is also true for any ab 
transitions: nn, neutrino oscillations. (The generalization for the relativistic case is 

trivial.) 

5.3 Cancellation of divergences in the potential model 

The consideration given above is infrared free. Let us return to the 5'-matrix problem for- 
mulation (oo, — oo). Due to the zero momentum transfer in the e- vertex, any matrix element 
contains the singular propagator (see Figs. 7 and 8a). However, the matrix element of potential 
model Tii obtained by means of 5*- matrix approach is not singular (see Eq. (9)). The same 
is true for the process shown in Fig. 7. From microscopic theory standpoint the reason is as 
follows. 

In addition to the singular propagator the matrix elements mentioned above also contain 
the block which is a sum of the zero angle rescattering diagrams of n. As a result, the self- 
energy part — V appears. The corresponding mechanism of the cancellation of divergences 
(the forming of the self-energy part) is illustrated by Eq. (19), where Gl^ ~ 1/0. 

We are interesting in off-diagonal matrix elements which do not contain the sum mentioned 
above (TJ^j instead of T^) and hence diverges because one singular propagator after e-vertex 
appears in any case. (Recall that the formal sum of series in e gives the meaningless self-energy 
part E ~ eVO.) 

The principal result of this section is that the FTA has been verified by the example of the 
exactly solvable potential model. It is involved in the block iT^{r) —<npO \ T^{t) \ Onp> as a 
special case. 

6 nn transitions followed by annihilation 

As shown above, the FTA reproduces all the potential model results. Besides, for non-singular 
diagrams it converts to the 5- matrix theory (see Sect. 6.1). We now proceed to the main 
calculation. 

Let us consider the process (2) in nuclear matter (see Fig. 8a). The Hamiltonians Hq and 
Hi{t) are the same as in Sect. 4. The 4-momenta of n and n coincide. The T'^-operator involves 
all the n-medium interactions. In consequence of this E = 0. In essence, we deal with 2-step 
nuclear decay: dynamical nn conversion, annihilation. Its dynamical part last only 10~^^ s. 
The sole distinction with respect to the decay theory is that the FTA should be used because 
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the antineutron propagator is singular. 
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Figure 8: (a) nn transition in the medium followed by annihilation, (b) Same as (a) but with 
escaping of particle in the nn transition vertex. 

We give the expressions for the amphtudes from Ref. [13]. Thereupon they will be obtained 
as a special case of a more general problem. The matrix element of the process shown in Fig. 
8a is 

Tf,{t) = -e 1^ dt,irf,{t,Q, (45) 

imt,Q = imr) =</| T^t) \0n,> . (46) 

Here TJ-- (r) is the matrix element of the antineutron annihilation in a time r — t — tc (compare 
with the matrix element of ^'-matrix (17)). The T'"(T)-operator is given by (39). Similarly to 
(40), the n-medium annihilation is separated out in the block 7Jj(T). 

Consider now the more general problem. We calculate the matrix element Tfi(t) shown in 
Fig. 8b on the interval (t/2, —t/2). As a result, it will be shown that: (a) If g 7^ (g is the 
4- momentum of particle escaped in the nn transition vertex) and i — > 00, we come to the usual 
^'-matrix amplitude, (b) If g — >^ 0, Eq. (45) is obtained. Such scheme allows to verify and 
study the FTA. Also we will see the point in which the standard calculation scheme should be 
changed. 

Consider the imaginary free-space decay 

n^n + ^, (47) 

$(x) = N^exp{—iqx), — {2qoQ)~^/^. For decay to be permissible in vacuum put = 
m — 2m$. As with Hnn, the decay Hamiltonian Ti'nn is taken in the scalar form T-i^n — 
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The corresponding process in nuclear matter is shown in Fig. 8b. This is a nearest analogy 
to the process under study. The neutron wave function is np{x) = Nnexp{—ipx), where — 
P = (po, p) ■ po = + p^/2m. The background nuclear matter field [/„ is omitted. 
Instead of Eq. (15) we have 

Hj = HU + H, (48) 
e' is dimensionless. In the lowest order in H'^^^ the matrix element Tfi{t) is 



Tfiit) = -< $J0 I Y.Tk{t) / ^ dt.KM |Onp>, 

dh... / dtkH(h)...H(tk). (49) 

-t/2 J-t/2 



k=i " 

rt/2 rt,,_ 

dh... 

-t/2 J-t/2 

Here < / | represents the annihilation products with (n) mesons. For the 3-tail 7Y^^ the relation 
Hnn \ Onp>= e \ Onp> used in Sect. 5.1, is invalid. The direct calculation is needed. 
Using the standard rules of quantum field theory, we obtain (see Appendix B) 

r*fc-i 

I -t/2 



Tfi{t) = -ie'NnN^ < /O I ^(-z)Tfc_i(t) / ~ dtj d''xkn\xk)e'^'^-'^>'^I{tk) \ 0>, (50) 

fe=i ■^-*/2 J 



^ / dt,- ^^^^-^kotk^M,o-Po+ko)_ (51^ 

-oo ZTT J-t/2 Kq — rUn — (p — q)^/2m^j + l\) 

From this point the calculations for Figs. 8a and 8b are essentially different. In Eq. (51) we 

put tk — oo and —t/2 — — oo. Then 

^ — / diee'*=(^°-^°+'=°) (52) 

27r ko-mn-{p- q)^/2mn + iOJ-oc ^ ^ 

and correspondingly 

7(oo) = G'e-^(P°-^°)*^ (53) 

^ Po-qo-mn-{p- q)2/2mn + iO ' ^^"^^ 
where G is the non-relativistic antineutron propagator. 
Let q — (0, 0) and irin — m (see Fig. 8a). Now 

G = ^— (55) 

Po — m — p^ /2m 

and ~ 1/0. This is an unremovable peculiarity. We deal with infrared divergence, what is 
obvious from Fig. 8a. We thus see the specific point (the limiting transition i ^ oo in (50)) in 
which the standard ^-matrix scheme should be changed. 
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6.1 Non-singular diagram 

Obtain now the amplitudes corresponding to Figs. 8a and 8b starting from (50). If g 7^ 0, the 
hmit t ^ 00 can be considered. In Eq. (50) we put t — > 00 and substitute Eq. (53). Taking 
into account that 

7V^e-*(^-^)^'= I >= ^n(xfe) I > (56) 
and using the relation / d^XkH\xk)'^n{,Xk) = H{tk) (see Appendix B), one obtains 

Tfi = -iN^e'G < /O I T" I Onp_g> . (57) 

Here | Ofip-q > is the state of the medium containing the n with the 4-momentum p — q, the 
T^-operator is given by (17). With the help of the relation 

iTf, =</| T |i>= N{27r)^S'ipf-pi)Mf, 

we rewrite (57) in terms of the amplitudes 

Msb = e'GM^. (58) 

Here Msb is the amplitude of the process shown in Fig. 8b, is the annihilation amphtude 
of n with the 4-momentum p — q, G is given by (54). We have obtained the usual ^'-matrix 
amplitude, which is the verification of (50). As in (16), the antineutron propagator is bare. 
It is easy to estimate the widths corresponding to Fig. 8b and free-space decay (47) : 

rs,«6'2r/(27r2), 

Ffree ~ e'2m$/(27r), (59) 
where we have put m$/m <^ 1. The ^-dependence is determined by the exponential decay law 

W^86,free = 1 - e-^^^-'f-^* - Fgft.freet- (60) 

These formulas will be needed below. 

6.2 Singular diagram 

Let g = and rUn — m (see Fig. 8a). In (50) one should put e' — e and — 1. Upon 
integration with respect to ko, Eq. (51) becomes 

I{tk) = -ie-^*'=^° /*' dtc. (61) 

J-t/2 
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As in (56), Nnexp{—ipxk) |0 >= ^fj(xfe) \ np>. Turning back to the Hamiltonian H{tk), one 
obtains 

°° i-tk 

Tfi{t) ^-ie<fO\J2 Tk{t) / dU I Orlp > . (62) 

Using the formula 

/ dti... dtk dtj{ti,...,tc)^ dtc dti... dtkf{ti,...,tc), (63) 

J-t/2 J-t/2 J-t/2 J-t/2 Jtc Jtc 

we change the integration order and pass on to the interval {t, 0). Finally 

Tfi{t) = -e /* dtc <fO\T^{t- U) I Onp >, (64) 

<mr\T) \Qn,>=iT%{T), 

which coincides with (45). The result is expressed through the submatrix TJ^{t). (Compare 
with (40).) Note that Tfi{t) coincides with the second term of (40) with the replacement 
<i |=<ripO|^</|. This can be considered as a test for the Tfi{t). 

Comparing Eqs. (64) and (57), one can see the formal correspondence: if g — > 0, GT^j — > 
ij^drmr). 

7 Infrared singularities and the formulation of the S- 
matrix problem 

In this section we consider the time-dependence of matrix elements and other characteristic 
features of the FTA and complete the calculation of process (2) (see also [13]). 

The FTA is infrared-free. It naturally connected with the conditions of experiment. Indeed, 
measurement of any process corresponds to some interval r. So it is necessary to calculate 
Ufi{T). The replacement U{t) — > S{oo) is justified if the main contribution gives some region 
A < T, so that Ufi{T > A.) — Ufi{oo) — Sfi =const. The expressions of this type are the basis 
for all ^'-matrix calculations. The following cases are possible. 

1. There is bound to be asymptotic regime. Then the usual scheme reahzed in field theory 
or non-stationary theory of scattering takes place. Fig. 8b corresponds to this case. 

2. There is no asymptotic regime. An example is provided by oscillation hamiltonian Hnn- 
We have usual non-stationary problem. The -S-matrix approach is inapplicable. Because of 
this, for Fig. 8a the calculation has been done in the framework of FTA. 

A somewhat different explanation of application of the FTA is as follows. If Tii = 7Y„n, 
the solution is periodic. It is obtained by means of non-stationary equations of motion and not 
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5'-matrix theory. This is clear from the ^'-matrix definition. To reproduce the hmiting case 
TC 0, i.e. the periodic solution, we have to use the FTA. 

Let us return to Eq. (64). The annihilation of n in nuclear matter can be considered as 
the decay of a one-particle state with the characteristic time Ta- Correspondingly, T^^ can be 
interpreted as the decay matrix of the n-medium state. Obviously 

77. (r > Ta) « Tf. = const. (65) 

and 

where VF" is the total decay probabihty of the n-nucleus. Let 

i > A fa Ta. (67) 

In view of this condition the submatrix Tj*^ can be calculated by means of S'-matrix theory. 
The FTA is needed only for description of the subprocess of the nn conversion. However, the 
condition (66) greatly simplifies the calculation. One can write immediately [13] 

Wait) ^ I -^^tTf^ r= eH^W = eH^ = Wf{t), (68) 

where Wa{t) is the probability of process (2). 

For nn transitions in nuclear Wt{t) — Wa{t) since all the n annihilate. The interpretation 
of Wa{t) has been given above: momentary nn conversion at some point in time between and 
t; annihilation in a time ~ 6 x 10~^^ s. The explanation of the t^-dependencc is simple. The 
process shown in Fig. 8a represents two consecutive subproccsscs. The speed and probability 
of the whole process arc defined by those of slower subprocess. Since ^ t, the annihilation 
can be considered instantaneous: for any ti < t the annihilation probability is W^"(t — ti) ~ 1. 
So, the probability of process (2) is defined by the speed of nn transition: Wa ~ Wf ~ t^, 
but not ~ t/r (see Eq. (6)). In essence, we deal with the limiting case r/t — > or, similarly, 
Tjj(r) = Tj^ at any r. Formally, the quadratic time-dependence follows from (64). 

Owing to annihilation channel, Wa is practically equal to the free-space nn transition prob- 
ability. So Tnn ~ Tnn, whcrc Tnn is thc osciUatiou time of neutron bound in a nucleus. 

All the results have been obtained by means of formal expansions. They are valid at any 
finite t. Consequently, the singularities of the ^S'-matrix amplitudes M^, and M2 result from the 
erroneus problem formulation. The problem should be formulated on the finite interval {t,0). 
If i — > 00, Eq. (68) diverges just as the modulus (16) squared does. The infrared singularities 
point to the fact that there is no asymptotic regime. 
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8 Summary and discussion 



The importance of unitarity condition is well known [19-21]. Nevertheless, the non-hermitian 
models are frequently used because on the one hand, they greatly simplify the calculation and 
on the other hand, it is hoped that an error may be inessential. This paper demonstrates that 
the non-unitarity of .S-matrix can produce a qualitative error in the results. Certainly, the 
unitarity is a necessary and not sufficient condition. We compare our results and potential 
model one. 

The time-dependence is a more important characteristic of any process. It is common 
knowledge that t-dcpcndcnce of decay probability in the vacuum and medium is identical. 
Equations (60) illustrate this fact. In our calculation the t-dependencies coincide as well: 
Wa Wt and Wf ~ t^. The potential model gives W^f"* ~ t, whereas Wf ~ t^. It is 
beyond reason to such fundamental change. 

The F-dependence of the results differs fundamentally as well. The probability of the decay 
shown in Fig. 8b is linear in F 

W^86 = ^mt ~ Vt (69) 
(see (59) and (60)). For Fig. 8a the annihilation effect acts in the same direction 

Wa^W F. (70) 

In the potential model the effect of absorption acts in the opposite direction Wf°^ ~ 1/F. 
Recall that the annihilation is the basic effect determining the process speed (see (6) and (68)). 
Let us consider the suppression factor R. From Eq. (68) we have 

Wn 

For similar processes the value ~ 1 is typical. Indeed, in the medium the free-space decay 
(47) suppressed by the factor 

^ (72) 



Ffree 7rm$ TT 

where we have put m$ ~ F. 

The realistic example is the pion production pn — > ppn' in vacuum and on neutron bound 
in a nucleus. If the pion energy is in the region of resonance, the pion absorption is very 
strong. This effects on the number of pions emitted from the nucleus, but not on the fact of 
pion formation inside the nucleus. (In the latter case the pion and products of pion absorption 
should be detected) . 

In the processes cited above it! ~ 1. The potential model gives it!pot — > 0: if F ~ 100 MeV 
and t ~ 1 yr [22], R^^t ~ 10~^°. 
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Consequently, in the potential model the t- and F-dependencies are principally incorrect. 
As a result, the suppression is enormous: i?pot — 0. This is not surprising since the potential 
model describes only Wn. Recall that in the strong absorption region Wn <^ Wa- 

The next important point is the competition between scattering and annihilation in the 
intermediate state. The models (13) and (15) reproduce the values of r and ri (see Sect. 4.3). 
Since the potential model does not describe the processes (2) and (11), it makes no sense to 
speak about competition effect in this model. The greater the |Imy|, the greater an error in 
the Wl°^ and Wa calculated by means of potential model. 

Consider now the effects of coherent forward scattering and absorption. Let there is a for- 
ward scattering alone: H — ReV. Since the FTA reproduces all the potential model results (see 
Sect. 5.2), it describes the above-mentioned special case as well, in particular, the suppression 
of oscillations by ReV. 

Let there is an annihilation vertex only: V — and 

n = Ha. (73) 

The annihilation Hamiltonian Tia is given by (27). In this case we inevitably arrive at the 
amplitude with singular propagator. The dressed propagator cannot arise in principal (see 
Sect. 4.2). In view of Eq. (22) the model (73) is reasonable and so the result Wa ~ Wf seems 
quite natural for us. In our calculation the approximation (73) has been not used. Nevertheless, 
the result is the same as in model (73). In this connection we briefly outline the principal points 
of our calculation. 

The process shown in Fig. 8b is described by the Hamiltonian Hj — H'^^ -\- H. Since H 
appears in the block T^- only, the antineutron propagator is bare. For Fig. 8a the picture is 
the same, however T/j ~ 1/0 (here we keep in mind the 5"- matrix problem formulation). Due 
to of this we had to use the FTA. 

The fact that antineutron propagator is bare is principal. It entails the divergence of the 
yS"- matrix amplitude; the application of FTA; the linear time-dependence of the matrix element 
Tfi{t) and i^-dependence of the result. In our opinion the models with dressed (and consequently 
non-singular) propagator are non- realistic (see Sect. 4). 

(Recall that in the potential model the antineutron propagator is dressed as E = y by the 
model construction. Since this model is inapplicable, the field-theoretical approach is used. 
The self-energy should be considered in the context of the concrete problem. Obviously, for 
Fig. 8b the propagator is bare. For Fig. 8a it is bare as well because the n-medium interaction 
is the same.) 

All the formulas up to (64) are true for any ah transitions in which = mj. (A generaliza- 
tions for the relativistic case and the case when nia 7^ m;, are simple.) The next important point 
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is the condition (67). For nn transitions in nuclei it is obvious because in this case the value 
t = To = 1.3 yr [22] is used (Tq is the observation time in a proton-decay type experiment). 
The condition t 3> leads to Eqs. (65) and (66). Due to of them the result does not depend 
on a specific form of H and coincides with the result given by the model (73). 

Once the antineutron annihilation amplitude is defined by (46), the rest of calculation 
is rather formal. The distinguishing features of the model is that the process amplitude is 
" propotional" to the annihilation amplitude T/j ~ T^j. This structure is typical for the direct 
processes. 

If the condition f ^ A is not fulfilled, the direct calculation of (64) is needed. However, the 
qualitative picture remains the same: the process amplitude is proportional to the absorption 
amplitude. 

It is interesting to study the behavior of Wt in the intermediate range t ~ A. It seems 
plausible that Wt depends shghtly on the value of A/t (in comparison with potential model 
results) . We also note that there is no asymptotic regime for free-space K^K^ oscillations. In 
our opinion, it makes sense to look at the calculation of Am = rriL — ms (GIM mechanism) 
from the standpoint of applicability of S'-matrix approach in this case (see Sect. 7). 

9 Conclusion 

The approach considered above reproduces all the results on the particle oscillations (Sect. 
5.2). Certainly, for the problems where the absorption is inessential, the standard model of 
oscillations is more handy since it is more simple. Our approach is oriented to the processes 
like (1) which are not described by the potential model. 

The direct calculation of nn transitions in nuclear matter followed by annihilation has been 
done. The results have been discussed in Sect. (8). We confirm our restriction [13] on the 
free-space nn oscillation time Tnn > 10^^ yr. Compared to [13], the result (68) was obtained as 
a special case of a more general problem. Besides, the medium corrections, the uncertainties 
related to amplitudes and competition between scattering and annihilation in the intermediate 
state have been studied. The model (73) and analysis made in Sect. 4 show that E = 0. 
Nevertheless, this is a point of great nicety. The further investigations are desirable. The 
region t < A and oscillations of another particles can be considered as well. 

The calculation up to (64) is formal. With the replacement T^i{r) — > Tfj^{r), where T^j 
is the 6-particle absorption amplitude, the matrix element (64) describes the process (1) in 
which nia — m;,. In this connection we point out some features of Eq. (64). a) The amplitude 
Tfi{t) is "proportional" to the amplitude Tj-{t). In the potential model the effect of 6-particle 
absorption acts in the opposite direction, which tends to suppress the process, b) In the lowest 
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order in e the potential model gives the linear t-dependcncc W^"^ ~ t/T. For any block T^^^r) 
model the time-dependence of the value | Tfi{t) p cannot be linear. 

Also we would like to emphasize that for a processes with zero momentum transfer the 
problem should be formulated on the finite time interval. 

Appendix A 

In this appendix the relative annihilation probability of the intermediate n for nn transition in 
the medium is calculated. Similarly to (31), we obtain the probability of process (2) in a unit 
of time 

= A^i 1 d$ I M2 p= a^N^ I I a^F, (74) 

flm = eC™. The normalization multiplier A^i is the same for r2 and F. The term "width" is 
unused because the t-dependence of process (2) does not need to be exp(— F2t) (see Sect. 7). 
In the low-density approximation [23,24] F = vpaa and 

F2 = al^vpaa. (75) 

The on-diagonal matrix element M™(0) corresponding to the process (n — medium) — > (n — 
medium) {n — medium) is 

M"*(0) = eC^Mf (O)G'^e = a^Mf (0) (76) 

(compare with (34)). Here M™(0) is the amplitude of zero angle scattering of n in the medium. 
Taking into account that 

;i2ImM-(0) = Ft, (77) 
-to 

l2ImM-(0) = vpat 
-to 

(To is the normalization time. Ft is the probability of the process (37) in a unit of time), one 
obtains 

Ft = al^vpat (78) 

and correspondingly 

^ = ri. (79) 
J- 1 

Equations (75) and (78) are interpreted in line with the low-density approximation physics. 
Appendix B 

The calculation is standard [25,26] up to the integration over t. The neutron and antineutron 
are assumed spinless. We have 

^n{x) \np>= ^n{x)b+{p) |0>= N^e-'P^ |0 >, (80) 
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< % I =< I AT^e^'^^ (81) 

Then 

< % I ^*{xc)'^n{xc) I np>^ 7V„7V^e^(^-PK. (32) 
In the last multipher of Eq. (49) we separate out the antineutron field operator "i/n^Xk)- 

H{tk) = J d'xkHixk) = J d^Xkn'{xk)^n{xk). (83) 

Equation (49) becomes 

00 /"tfc— 1 r 

Tfi{t) = -ie'NnN^ < fO I J2i-i)Tk-i{t) dtk j d'xkH' {xk)J{tk) \ >, (84) 



J{tk) = 1'^^ dU j d^x, < T{^n{xk)^n{xc)) > e'^^-^'^^^. (85) 

For Fig. 8a the problem is non-relativistic and so for Fig. 8b we also take the non-relativistic 
antineutron propagator 



< T{^lJn{xk)%{xc)) >= tG{xk -x,)=i j ^e-'=°(*'=-*=) j — - 



(27r)3 ko-mfi- \i^/2mn + iO 

(86) 

Upon integrating over Xc and k we obtain Eqs. (50) and (51). 
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